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Abstract 

The scalar Yukawa model, in which a complex scalar field, tp, interact via a real scalar field, %, is reduced 
by using covariant Green functions. It is shown that exact few-particle eigenstates of the truncated QFT 
Hamiltonian can be obtained in the Feshbach-Villars formulation if an unorthodox "empty" vacuum state is 
used. Analytic solutions for the two-body case are obtained for massless chion exchange in 3+1 dimensions 
and for massive chion exchange in 1+1 dimensions. Comparison is made to ladder Bethe-Salpeter, Feymnan- 
Schwinger and quasipotential results for massive chion exchange in 3+1. Equations for the three-body case 
are also obtained. 
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1. Introduction 



The model scalar field theory, based on the Lagrangian density (h = c = 1) 

C = d u tp* (x) d u (p(x) — m 2 Lp* (x)ip(x) 

+ \d v x{x)d v \{x) - ^ 2 x 2 (x)- g<P*{x)<p(x) X {x) (1) 

is often used as a prototype QFT in many studies. In particular, for the case [i = 0, it is known as the Wick- 
Cutkosky model [1,2]. This model has been investigated in various formalisms, in addition to the original 
ladder Bethe-Salpeter formulation [3], including the light-cone formulation [4,5], and others (see refs. [6-9] 
and [18, 20-23] and citations therein). 

We shall consider a reduced form of this theory in the present paper, in which the mediating chion field is 
partially eliminated by means of covariant Green functions, in a manner analogous to that discussed recently 
for QED [10,11]. In addition, we shall use the Feshbach-Villars (FV) formulation [12] for the complex phion 
field, and an unconventional "empty" vacuum state, as has been used recently to determine solutions of the 
\(<p*(p) 2 theory [13]. 

The fields ip and x °f tne model (1) satisfy the equations 

d v d vX {x) + ti 2 x(x) = -g<p*(x)<p{x), (2) 
d v d v tp{x) + m 2 ip(x) = -gip(x)x(x), (3) 

and its conjugate. 

Equation (2) has the formal solution 

X(x) = X o(x) + J dx' D(x - x 1 ) p(x'), (4) 

where dx = d N xdt in N + 1 dimensions, p(x) — —gip*(x)(p(x), Xo(x) satisfies the homogeneous (or free 
field) equation (eq. (2) with g — 0), while D(x — x') is a covariant Green function (or chion propagator, in 
QFTheoretic language), such that 

{d u d v + n 2 )D{x-x') =S N+1 {x-x'). (5) 

Equation (5) does not specify D(x — x') uniquely since, for example, any solution of the homogeneous 
equation can be added to it without invalidating (5). Boundary conditions based on physical considerations 
are used to pin down the form of D. 
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Substitution of the formal solution (4) into cq. (3) yields the equation 

8»8 M x) + mMz) = -M*)XoW - Mx) J dx'D(x - X 'W). (6) 

Equation (6) is derivable from the action principle 5 j dx C = 0, corresponding to the Lagrangian density 

C =d v ip* (x)d„<p(x) - m 2 tp* (x)<p>(x) - gp* {x)ip(x) xo {x) (7) 
+ i J dx'p{x)D{x - x')p(x'), 

provided that D(x — x') = D(x' — x). 

The QFTs based on (1) and (7) are equivalent in that, in conventional covariant perturbation theory, 
they lead to the same invariant matrix elements in various order of perturbation theory. The difference 
is that, in the formulation based on (7), the interaction term that contains the propagator corresponds to 
Feynman diagrams involving virtual chions, while the term that contains xo correspond to diagrams that 
cannot be generated using the term with D(x — x'), such as those with external (physical) chion lines. 

We shall consider this scalar theory in the Feshbach-Villars (FV) formulation [12]. The reason for 
doing so is that this leads to a QFThcorctic Hamiltonian which is Schrodingcr-likc in form, for which exact 
eigensolutions can be readily written down. In the FV formulation, the field ip and its time-derivative ip are 
replaced by a two-component vector 



u = —r= (nvp + i<p) 
(mtp — ip) 



v = 



/2m 



so that, for example, 2mip*<p — (u* + v*)(u + v) = (^r\T<\>, where r\ and r are the matrices 

and t 



1 

-1 



1 1 

-1 -1 



In the FV formulation the equation of motion (3) takes on the form 



or, upon using (4), the form 



^ = _ i V2r(/) + m ^ + 2^ r0X ' 



2 ^"i VV(/) + m ^ + 2^ rto 



(8) 



(9). 



(10) 



+ — T(j) / dx' D(x - x')p(x'), 
2m J 

where p = —gtp*tp = ———<$ rjT(j). Equation (11) is derivable from the Lagrangian density 



(11) 



2m 1 



Cfv{ x ) = i<t>^ {x)i]<j>(x) — - — V<p(x) ■ \7(f>(x) — m(fJ (x)(j)(x) 

2m 



p(x)D(x — x')p(x'), 



(12) 



where = (f^rjr. Note that £ of eq. (7) is not identical to Cfv- Indeed C = Cfv + Jj(<£*</i>). However, they 
lead to identical equations of motion ((6) and (11)), and so are equivalent in this sense. Henceforth, we base 
our results on Cfv- 

We note that the momenta corresponding to u and v are 

ac . . 

Pu = ttt = iu , and p v = -iv , 
ou 

that is, u* and v* are, in essence, the conjugate momenta, so that the Hamiltonian density is given by the 
expression 

H(x) = ^(x) V h(x)(f>(x) + £^(x)<b(x) X o(x) -\j dx'p(x)D(x - x')p(x'), (13) 

where h(x) = r (^2m)V 2 +mr]. 

We use canonical equal-time quantization, whereupon the non-vanishing commutation relations are 

[u(x,t),p u (y,t)] =i6 N (x-y) and [«(x, t),p v (y, t)} - i6 N (x - y), (14) 

or equivalently, 

[0a(x,t),4(y»*)]='?a6<5 JV (x-y), a, 6 - 1, 2 (15) 

and where T = [<f>i — u, (f>2 — v], while r\ a h are elements of the r\ matrix (9). Using these commutation 
relations, the QFTheoretic Hamiltonian can be written as 

H = J d N x[Ho(x)+H x (x)+H I {x)], (16) 



where 



Ho(x) =^(x)r 1 h(x)<j)(x), (17) 
H x (x) = JL-$( x )cP(x)xo(x), (18) 



and 



H i( x ) = -J^ I dx' (t>{x)4>{x)D{x - x')(j){x')4>{x') 



8m 2 

.,2 



J dx' D(x-x')(l>(x){(f>{x')(f>(x'))(t>(x), (19) 



8m 2 

and where we have used r 2 = in the last step of (19). Note that no infinities are dropped upon normal 
ordering, since none arise on account of the r 2 = property. 

As already mentioned, Hi contains the covariant chion propagator, hence in conventional covariant per- 
turbation theory it leads to Feynman diagrams with internal chion lines. On the other hand, H x corresponds 
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to Feynman diagrams with external chions. However, we shall not pursue covariant perturbation theory in 
this work, and so shall not consider that approach further. Rather, we shall consider an approach that 
leads to some exact eigenstates of the Hamiltonian (16), but with H x = 0. Such a truncated Hamiltonian 
is appropriate for describing states for which there is no annihilation or decay involving the emission or 
absorbtion of real chions. 

In the Schrodinger picture we can take t = 0. Therefore, we shall use the notation that, say </>( x , t = 
0) = </>(x), etc., for QFT operators. This allows us to express (19) as 

W 7 (x) = J d Ar S 'G(x-x')^(x)(^x')0(x'))^(x), (20) 

where 

G( "WI D( " V *' = csf/ """""^'wtj?- ,21) 

Explicitly, for N = 3 this becomes 

1 p -Ai|x-x'| 



for N = 2 it is 



G(x-x') = ^o(M|x-x'|), (23) 



where Kq(z) is the modified Bessel function, whereas for TV = 1 it has the form 

G(x - x') = -U-H*-*'! . (24) 
At this stage we proceed in an unorthodox fashion, and define an empty vacuum state, |0), such that 

0a|6)=O. (25) 

This is different from the conventional Dirac vacuum |0) (the "filled negative energy sea" vacuum), which 
is annihilated by only the positive frequency part of ip and by the negative frequency part of ip* (see, for 
example, ref. [14], p. 38). With the definition (25), the state defined as 

|Vi>= |^t(x>7/(x)|0), (26) 

where f(x) is a two-component vector, is an eigenstate of the truncated QFT Hamiltonian (16) (H x = 0) 
with eigenvalue E\ provided that the /(x) is a solution of the equation 

/i(x)/(x) - ^/(x). (27) 

This is just the free-particle Klein-Gordon equation for stationary states {\ipi} is insensitive to H\), in the FV 
formulation . It has, of course, all the usual negative-energy "pathologies" of the KG equation. The presence 



of negative-energy solutions is a consequence of the use of vacuum (25). However, that is the price that has 
to be paid in order to obtain exact eigenstates of the truncated Hamiltonian ( eq. (16) with H x — 0). We 
shall refer to as a one-KG-particle state. 



2. Two-particle eigenstates 

We can define two-KG-particle states, analogously to (26): 



- / d N xd N y F a6 (x,y) 4(x)4(y)|6>, (28) 



where summation on repeated indices a and b is implied. This state is an eigenstate of the truncated QFT 
Hamiltonian ( (16) with H x = 0) provided that the 2x2 coefficient matrix F = [F a b] is a solution of the 
two-body equation, 

r ? /i(x)r ? F(x, y) + [ V h(y) V F T (x, y)] T + V(x - y)r T F(x, y)r = £ 2 F(x, y), (29) 

where the superscript T stands for "transpose" . The potential here is given by 

nx-y) = -^G(x-y), (30) 

where G is specified in equations (21)-(24). Equation (29) is a relativistic two-body Klein-Gordon-Feshbach- 
Villars-like equation, with an attractive Yukawa interparticle interaction. If V = 0, then cq. (29) has the 
solution F(x, y) = 5i(x)^J(y), where each /»(x) = r?<7i(x) is a solution of the free KG equation (27), with 
eigenenergy £j, where E2 = £1 + £2, as would be expected. 

In the rest frame, P total 1^2) = 0, equation (29) simplifies to 

h(r)F(r) + [h(r)F T (r)] T + V{y)t t F(y)t = E 2 F(r), (31) 

g 2 

where r = x y, h = ntvn, and VYr) = — — — rGfr) in this case. It is useful to write out this equation in 
component form, with 

s(r) t(r) 



F(v) 
namely 



u(r) v(r) 



(32) 



--3-V 2 (2s -u-t) + V(s-t-u + v) + (2m - E 2 )s = 0, (33) 

Zi 77?- 

--?-X7 2 (s-v) + V(s-t-u + v)-E 2 t = 0, (34) 
--^-V 2 (s -v) + V(s - t - u + v) - E 2 u = 0, (35) 

Zi 77?- 



and 



2m 



V 2 (t + u - 2v) + V(s-t-u + v)- (2m + E 2 )v = 0. 



Equations (34) and (35) imply that t(r) = u(r), so that only three equations survive: 

(2m -E 2 - — V 2 + V)s + ( — V 2 - 2V)t + Vv = 0, 
m m 

(-^ + V) S -(E 2 + 2V)t + (±Vi + V) V = 0, 



and 



Vs + (-—V 2 - 2V)t - (2m + E 2 - —V 2 - V)v = 0. 
m m 



(36) 

(37) 
(38) 

(39) 



These equations have positive-energy solutions of the type E 2 = m + m + ■ ■ ■, negative-energy solutions of 

the type E 2 = — m — m H , and "mixed" type solutions with E 2 = m — m-\ (this is clear, for example, 

if V — and the particles are at rest). 

For the positive-energy solutions, if we write E 2 = 2m + e, then in the non-relativistic limit |(e + V + 

p 2 

— )v\ <C \mv\ (and similarly for s and t), and so equation (38) and (39) show that t and v are small and 
m 

doubly-small components respectively, by factors O ) • Thereupon, equation (37) reduces to 



--V 2 .s(r) + V»s(r) =es(r), 
m 



(40) 



which is the usual time-independent Schrodinger equation for the relative motion of two particles, each 
of mass m, interacting through the potential V(r). Similarly, in the non-relativistic limit, v is the large 
component for the negative-energy solutions ( i.e. E 2 = —(2m + e), and s — > v, V — > — V in (40)), while 
t is the large component for the mixed energy solutions. This is obvious from the form of the free-particle 
solutions (V = 0), which are 



F(r) = s 
for E 2 = 2lu = 2 v / p 2 Tm 2 , 



for E 2 = 0, and 



F(r) = v Q 



i {^-Y 

/ V \2 ( uj — m \2 



F(r) = t 



(^) 2 i 

v oj-\-m ' 



2- < 1 



so 



1 



( _P_ \2 ( P \4 

\2m> V 2m / 



P 

2m 2 +p 2 
1 



2m 2 +p 2 _ 



< 1 



^0 



V 2m / V 2m > 



o l P r 



(41) 



(42) 



(43) 



for i? 2 = — 2uj = —2y / p 2 + m 2 , and where so,io and wo are constants. 

Equations (37)-(39) can be reduced by taking suitable linear combinations, whereupon it follows that 
(E = E 2 ) 

(2m + E)v = (2m - E)a + 2Et (44) 
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and 

[E(4m 2 - E 2 ) - 8m 2 V]s = -[E(2m + E) 2 + 8m 2 V]t. (45) 

It is easily verified that the free particle solutions (41) - (43), in particular, satisfy these relations. One can 
therefore write 

s = [8m 2 V + E{2m + E) 2 ]w (45) 

and 

t = [8m 2 V + E(E 2 -4m 2 )]w, (46) 

where w is a solution of (E ^ 0) 

--V 2 w+-^-Vw = -((E/2) 2 -m 2 )w. (47) 
m E/2 m 

Once (47) is solved for w, the components s,t = u, v of the matrix F follow from (45) and (46). 

Equation (47) is form-identical to the Schrodingcr equation, and so can be solved in the same manner as 
the latter for both bound and continuum states. In general this has to be done numerically. In some cases, 
such as for fj, = (massless chion exchange) in 3+1, and for fi ^ in 1+1, analytic solutions of Schrodinger's 
equation are known. We shall not discuss solution of equation (47) for the entire range of the parameters \x 
and g, for various N. Rather, we shall consider bound states for the two analytically solvable cases in some 
detail, leaving most of the rest for another time. We shall also examine bound states for one 3+1 case, with 
n/m = 0.15, numerically, since this case was studied in some detail recently by Nieuwenhuis and Tjon [8b]. 



3. Two-body bound states in 3+1 for massless chion exchange. 

We consider, first, the solution of equation (47) for N — 3 and massless chion exchange (i.e. fi = 0). In 
this case one can use the known hydrogenic solutions of the Schrodingcr equation to obtain the solutions of 
cq. (47). Thus, for the bound states we obtain the result 

(W-m 2 = -^, (48) 

g 2 

where a = » , and n — 1, 2, 3, • • • is the principal quantum number. This yields the positive energy 

two-particle bound-state spectrum 

*- J 2 ( 1+ ^(3>™H© 2 -A© 4+ -) • m 



which is seen to have the correct Rydberg non-relativistic limit. Note that the relativistic spectrum retains 

the "accidental" Coulomb degeneracy with respect to £, unlike the Klein-Gordon spectrum for a static 

electromagnetic potential. This may seem surprising at first glance, but upon reflection it is not so. The 

potential V = —eA° = — — enters both linearly and quadratically for the KG equation in an external 

r 

electromagnetic field = (A° 7 0), whereas in the present scalar Yukawa theory the potential enters only 
linearly . 

The two-particle bound-state wave functions corresponding to the eigenenergies (49) can be lifted simi- 
larly from the Schrodingcr hydrogenic results. For example, the ground state wave function (unnormalized) 

corresponding to (49) with n = 1 is w = e~^ r ', where j3 — ——. The critical value of a beyond which E (eq. 

E 

(49)) ceases to be real is, evidently, a = 1 (a = n, in general). 

The question arises how these positive-energy bound-state solutions compare to corresponding results 
obtained in other formulations of this model. The original Wick-Cutkosky solutions [1,2] of the massless 
chion case in the ladder Bethe-Salpeter approximation, as well as the corresponding light-cone formulation 
[4,5], give the small-a expansion 



This is different from the present results (49), for which the lowest-order correction to the Rydberg energy 
is 0(a 4 ). The unusual a 3 In a and a 3 terms (which have been termed a "disease" of the ladder Bethe- 
Salpeter Wick-Cutkosky solution [15]) are, apparently, an artifact of the ladder approximation. Variational- 
perturbative calculations of the scalar Yukawa model [9] that use a conventional (Dirac "filled negative-energy 
sea") vacuum, also yield no 0(a 3 In a, a 3 ) terms. However the 0(a 4 ) terms of [9] are different from those of 
equation (49). For example, for the ground state [9] give 



and there is no ^-degeneracy for states with n > 1. The disagreement between (51) and (49) at 0(a 4 ) may 
be, in part, because virtual annihilation effects, which contribute at 0(a 4 ), have not been included in the 
calculations of ref. [9] (i.e. eq. (51)). We might mention that a comparison of corresponding results for the 
Coulomb QED model (i.e. QED in the Coulomb gauge with the transverse a ■ A interaction turned off), 
namely calculations analogous to the present that use an "empty" vacuum [16], and conventional- vacuum 
variational-perturbative results [17] are in agreement at 0(a 4 ). Thus, the reason for the lack of agreement 
between (49) and (51) beyond 0(a 2 ) in the perturbative domain (a <C 1) is, at present, not clear to the 
author and is a matter that needs further investigation. The disagreement with the ladder Bethe-Salpeter 




(50) 




(51) 
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results (50) beyond 0(a 2 ) is less surprising, given the oft-mentioned shortcomings of the ladder Bcthc 
Salpeter approximation (see, for example, the discussion in refs. [8b], [15], [21], [24]). 

It might be tempting to speculate that the reason for the disagreement in the perturbative domain lies 
buried in the use of the unconventional vacuum (25). However, as has already been mentioned, the use of 
such a vacuum in the Coulomb QED case, leads to no disagreement up to 0(a 4 ). The present approach 
(with an empty vacuum) seems to correspond to the summation of a particular subset of diagrams in the 
conventional perturbative treatment (namely ladder and crossed-ladder diagrams; see also sect. 5). Such 
has been shown to be the case for the model theory where a (second) quantized Dirac field interacts with a 
classical (c-number) electromagnetic field [30]. However, since a detailed analysis along the lines of [30] has 
not been carried out for the scalar Yukawa model, this remains a point of speculation at this stage. 

In the non-perturbative regime, the present results (which are not perturbative) decrease monotonically 

with increasing a to E = \/2m at — = 1. This behaviour is characteristic of calculations of E(a) for 

n 

this model in that E(a) decreases monotonically from E(a = 0) = 2m in all cases. However, the various 
approaches yield results that differ markedly in detail. For example, for the ground state (n = 1), the present 
results decrease much more rapidly with increasing a than any of the ladder Bethe-Salpeter calculations 
[2,4,5,6,8] or the Haag-expansion results of Raychaudhuri [18]. In particular, none of the latter give a 
restriction < a < 1 as does the present analytic formula (49), and so no critical value of E(a = 1) = \[2m. 
Indeed the ladder Bethe-Salpeter and the variational -perturbative approximations all give values ofE/m, for 
n=l, which are near or above 1.9 at a = 1, far above the present value of 1.414 • • •, though the variational- 
perturbative results fall below the ladder Bethe-Salpeter ones for strong coupling. (A comparative plot of 
the ladder Bethe-Salpeter and variational- perturbative values is given in fig. 1 of ref. [9]). Raychaudhuri's 
Haag-expansion results for fi = are also very different in detail from the present analytic results (see table 
II of ref. [18]), however they do exhibit a critical value of a, but at a much larger value of a^ ay ~ 2.31 
(versus a c = 1 in our case) with the corresponding value of E^ &y — 1.30 (vs. 1.414 • • • here). 

There are no negative-energy bound state solutions in the present system since, as has been mentioned 
earlier, the potential effectively reverses sign for the negative-energy case (this also happens in the case 
of one-particle relativistic equations, such as the Klein-Gordon-Coulomb and Dirac-Coulomb equations). 
However, eq. (48) has "mixed-energy" type solutions with 



These unphysical solutions do not have a Rydberg non-relativistic limit, and arise because of the retention 
of negative-energy solutions in the present formalism. For these mixed-energy solutions E increases mono- 



E = m 



(52) 
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tonically with increasing a from a value of E — at a = to the value E/m = \[2 at a = n. It is of interest 
to note that the positive-energy and mixed-energy solutions join smoothly at a = n. Thus, for < a < a c , 
E(a) forms a continuous double-valued function, with the upper branch being the positive-energy solution 
and the lower branch being the mixed-energy solution. This is clear if one notes that eq. (48) can be recast 
into the equation of the semicircle I - — ^ — 1 J + y— J = 1, a > 0. Thus, eqs. (49) and (52) correspond, 
respectively, to the upper and lower branches of this semicircle. Solutions that include such "mixed-energy" 
behaviour arise in some other formulations of the relativistic two-body system. In particular Raychaudhuri's 
Haag-expansion results exhibit this "double-valued" behaviour (see fig. 5 of ref. [18]), although he does not 
identify the lower (in energy) branch as a "mixed-energy" phenomenon. 



4. Two-body bound states in 1 + 1. 

Equation (47) in 1+1 corresponds to a one-dimensional Schrodingcr equation with an exponential po- 
tential. This happens to be one of the not numerous cases for which the bound state eigenvalues can be 
expressed in terms of common analytic functions [19]. Thus, for the present case, the bound state eigenvalues 
are given by 

JKl) = even parity, and J„(7) = odd parity, (53) 

where J v {"f) is the usual Bessel function, while v = 2\Jm 2 — (E/2) 2 //i and 7 = gj ' yj (j?E. We evaluated some 
solutions of equation (53) using the Maple programme. Their general behaviour is similar to the expressions 
(49) and (52), in that, for given m and \i there are positive-energy solutions, for which E/m starts from 2 
when g = and decreases monotonically with increasing g for < g < g c , as well as mixed-energy solutions 
for which E/m start from at g = and increases monotonically with increasing g for < g < g c . The 
two curves become coincident at E(g c ), as was the case for the 3+1 solutions (49) and (52). A sample of 
the positive-energy ground-state solutions of (53), for m/fi = 6.944, is given in table 1. We also list, for the 
ground state, corresponding non-relativistic results as well as perturbative and variational results obtained 
previously for the real scalar Yukawa model in 1+1 (i.e. real phions exchanging real chions) [19]. The results 
for E are listed in the table for various values of the parameter A = g / {A^/itm) . This parameter is not 
dimcnsionlcss in 1+1, unlike in 3+1 (where a = A 2 is dimcnsionless). The present relativistic results are 
seen to fall increasingly below the variational and perturbative results as A increases. This is reminiscent of 
what was observed in 3+1. Also, the variational results (like the perturbative and non-relativistic ones) give 
no indication of a critical value of A c = g c /(4:y/irm) beyond which no real solution for E is obtained from 
(53). 
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The behaviour of the first few even-parity excited-state solutions of eq. (53) is qualitatively similar 
to that of the ground state. Table 2 is a list of the critical values of E(\ c ) for the ground and first three 
even-parity excited states, for m/ fi = 6.944. 

The mixed-energy solutions of eq. (53) for the ground state of the mj /i — 6.944 case are listed in table 
3. We have not done an investigation of the spectrum of energies for the entire range of values of mj \i in 
1+1, as these can be readily obtained from eq. (53) as needed. 

5. Two-body bound states in 3+1 for massive chion exchange. 

Nieuwenhuis and Tjon recently reported an interesting study of the ip 2 x model using a Feynman - 
Schwinger formulation [8b] that contains all ladder and crossed ladder diagrams. They give a detailed com- 
parison of their results with ladder Bethe-Salpeter and various "quasipotential" equations (i.e. modifications 
of the Bethe-Salpeter equation), among them the Blankcnbecler-Sugar equation [20], the Gross equation [21] 
and the "equal-time equation" [22,23]. Since these various equations are written out and discussed in ref. 
[8b] we shall not repeat this here. Nieuwenhuis and Tjon find that the ladder Bethe-Salpeter results increas- 
ingly underestimate the two-body binding energy, as a increases, compared to their Feynman-Schwinger 
ladder-plus-crossed-ladder results, and, indeed, compared to the various quasipotential results (see fig. 1 of 
ref. [8b], in which E/m are plotted for < a < 0.93). The various quasipotential results are distributed 
between the ladder Bcthc-Salpctcr and the Nieuwenhuis and Tjon Feynman-Schwinger results. 

It is of interest, therefore, to compare the predictions of the present approach with those given in the work 
of Nieuwenhuis and Tjon [8b]. Thus, we calculated the positive-energy E/m values for various < a < a c 
for this n/m = 0.15 case. We used the Maple differential equation solver and the "shooting method" (see, 
for example, ref. [25] sect. 2.5). We tested the accuracy of this numerical procedure on the analytic cases 
discussed above in sections 3 and 4, and found that the analytic results could be reproduced with essentially 
arbitrary accuracy (we tried up to 10 figures for some cases). Since the results of the various approaches 
compared by Nieuwenhuis and Tjon , as well as ours, converge at small a (to the non-rclativistic values), 
and since they all exhibit the same monotonically decreasing shape, we present here a list of E/m only at 
a = 0.9 where the differences are quite marked. These results, in decreasing order of E/m, arc: Ladder 
Bethe-Salpeter 1.963, Blankenbccler-Sugar 1.931, Gross 1.910, present results of eq. (47) 1.882, Gross (with 
retardation) 1.880, equal-time 1.860 and Nieuwenhuis and Tjon 1.770. These values, except for the present 
calculation, were read off Fig. 1 of Nieuwenhuis and Tjon [8b] and so the quoted accuracy of the last figure 
is doubtful. A more detailed list of our results for this \i/m = 0.15 case is given in Table 4, alongside the 
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Nieuwenhuis and Tjon and Gross-equation results (taken from Fig.l of [8b]). 

A pictorial representation of most of these positive-energy solutions is given in Fig. 1. We plot our 
present results (solid line), along with the Ladder Bethe-Salpeter (dashed curve), the Blankenbecler-Sugar 
(diamonds), Gross with retardation (crosses), equal-time (stars) and Nieuwenhuis and Tjon (circles) values. 
All these results, save ours, are taken from Fig. 1 of Nieuwenhuis and Tjon [8b]. 

As can be seen, our results are most similar to those of Gross (with retardation). As such, they are 
substantially below the Ladder Bethe-Salpeter results, and close (or closer) to those with ladder and cross 
ladder diagrams. This suggests that the present formalism is equivalent, in some sense, to a summation of all 
ladder and crossed ladder diagrams in covariant perturbation theory. The fact that the present formalism is 
non-perturbative and contains the chion propagator as the effective potential, is supportive of this conjecture. 
However, it would be necessary to do an analysis like that of Guiasu and Koniuk [30] to demonstrate this 
explicitly and convincingly. 

It is interesting to note that the Gross results exhibit the "double- valued" structure of E{a), with an 
upper "positive-energy" and a lower "mixed-energy" branch that join continuously at a G (see inset in figure 
1 of Nieuwenhuis and Tjon [8b]). This is precisely the behaviour that arises in the present formalism (e.g. 
eqs. [49] and [52] for /i = 0). Indeed, the Gross critical value of E(a G = 1.28) ~ 1.4m is quite close to the 
value that we obtain, namely E(a c = 1.2087) = 1.48386m. Our mixed-energy results for fi/m = 0.15 are 
also similar to those obtained from the Gross eq., though our results lie somewhat higher. For example, we 
obtain E(a = 0.25) = 0.219, E(a = 0.5) = 0.446, E(a = 1.0) = 0.985, whereas the Gross-cq. results are 
E G (a = 0.25) = 0.07, E G (a = 0.5) = 0.36, E G (a = 1.0) = 0.83. 

Table 5 contains a list of our E(a) for fi/m = 0.15 in 3+1 for the two lowest excited states (labelled 
n = 2) with i = (radial wave-functions w(r) with one node) and 1=1 (nodeless w(r)). The general shape 
of E(a) for these excited states is similar to that of the ground state, except that the critical values of a 
(which are similar, but not coincident for these n = 2 states) are over twice as large as for the ground state. 
This is reminiscent of the a c = n behaviour of the massless chion case (p = 0). There are mixed-energy 
branches of E(a) for the excited states, just as for the ground state (and the analytic /i = case), however 
we do not give a list of them here. Their behaviour is very similar to the case already discussed, in that 
E(a) increases monotonically from E(a = 0) = to E(a c ). 
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6. Three-body eigenstates 



It is straightforward to write down three-body states analogous to the two-body state (28), namely 

|Va> = J d N x 1 d N x 2 d N x 3 F abc (yL 1 , X2 ,X3)(pi(xi)4(x2)4(x3)\U}- (54) 

These are eigenstates of the Hamiltonian (16) (with H x = 0) corresponding to the eigenvalue E 3 provided 
that the 2 3 = 8 coefficient functions F a f, c (xi, x 2 , x 3 ) are solutions of 

ft, afe (xi)F fe&c (xi,X2,X 3 ) + /j bfe (x 2 )F afec (xi,X2,X3) + /l cfe (x 3 )i ;1 af)fc (xi,X2,X3) 

+ V(xi - x 2 )f afcl f bfe2 F fclfc2C (xi,X2,x 3 ) + V(x 2 - x 3 )f bkl f ck2 F aklk2 (x 1 ,x 2 ,X3) (55) 

+ ^(x3-Xi)f cfel f a fe 2 F fel (, fc2 (xi,X2,X3) = £ 3 F a (, c (xi,X2,X 3 ) . 

Equations (55) are relativistic three-body Klein-Gordion-Feshbach-Villars-like equations. They are similar in 
structure to three-fermion equations derived previously [11, 26, 27], except that h is then the Dirac operator 
and a, 6, c, ... are Dirac spinor indices, so that there are 4 3 = 64 equations (not all of them independent, 
however). Solution of the equations (55) is much more challenging than for the two-body case, as it is for 
any relativistic three-body problem [28, 29, 31]. We shall not discuss the solution of these equations here. 
Generalizations for iV-body eigenstates can be written down in an analogous fashion. 

Concluding remarks 

We have shown that the scalar Yukawa model can be recast in a form such that exact few-body eigen- 
states of the QFTheoretic Hamiltonian, in the canonic equal-time formalism, can be determined for the 
case where there are no free (physical) quanta of the mediating field (i.e. only virtual quanta). This is 
achieved by the partial elimination of the mediating field by means of Green functions, as well as by the 
use of the Feshbach-Villars formulation of scalar FT and the use of an "empty" vacuum state. This last 
requirement leads to the retention of negative-energy solutions, akin to the one-particle Klein-Gordon and 
Dirac equations. 

We considered the solution of the resulting two-particle equations in some detail, for massless and 
massive quantum exchange, in 1+1 and 3+1 dimensions. Analytic solutions for the two-particle bound state 
eigenenergies were obtained for massive exchange in 1+1 (eq. 53), and for massless exchange in 3+1 (eq. 
49). For the massive exchange case in 3+1 our results compare favourably with recent covariant Bethe- 
Salpeter based models, which include ladder and crossed-ladder diagrams ( particularly those of the Gross 
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equation). However, our results give stronger two-body binding energies than results which contain only 
ladder diagrams (such as the ladder Bethe-Salpeter results), but weaker binding than results that go beyond 
ladder-plus-crossed-ladder effects (such as the Feynman- Schwinger calculations of Nieuwenhuis and Tjon ). 

The present approach can be used for relativistic three-body systems, and we derived such equations for 
the present scalar Yukawa model. It can also be used for other QFT models, such as spinor Yukawa model 
and QED. 
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Table 1. Values of E 2 (ground state) for in = 6.944 (units: /i = 1) in 1+1. 



A 


E 2 


E 2 


E 2 [19] 


E 2 [19] 


g 

4y / 7rm 


eq. (53) 


non-rel. 


variational 


pert. 


0.0 


13.888 


13.888 


13.888 


13.888 


0.05 


13.8866 8289 


13.8866 8318 


13.888 


13.887 


0.1 


13.8746 2814 


13.8746 5453 


13.879 


13.875 


0.2 


13.7885 3427 


13.7898 3995 


13.805 


13.791 


0.4 


13.2754 9678 


13.3205 9729 


13.384 


13.334 


0.6 


12.0275 4474 


12.4161 0158 


12.578 


12.472 


0.7 


10.7258 6117 


11.7945 8644 


12.578 


12.472 


0.75 


9.2844 6172 


11.4408 8286 


12.578 


12.472 


0.7589 10218* 


8.2861 8251 








0.8 




11.0583 6849 


11.404 


11.201 


1.0 




9.2380 4702 


9.887 


9.528 


1.2 




6.9492 3608 


8.057 


7.454 


1.4 




4.1878 1169 


5.940 


4.985 


1.6 




0.9507 0057 


3.561 


2.122 


1.8 




-2.7644 9139 


0.942 





Table 2. Values of E 2 (\ c ) in 1+1 for m = 6.944 (units: ^ = 1) for the ground state (labelled n=l) and first 
three even-parity excited states (labelled n=3,5,7). 



n 


A c 


£2 (A c ) 


E 2 (X c )/m 


1 


0.7589 1021 83 


8.2862 


1.1933 


3 


1.0429 3738 8 


8.9174 


1.2842 


5 


1.2762 9938 2 


9.3428 


1.3454 


7 


1.4983 1853 8 


9.6932 


1.4348 
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Table 3. Ground-state mixed- 



energy solutions E 2 (X) in 1+1 for m = 6.944 (units: \i = 1). 



A 


E 2 {\) 


0.1018 5940 83 


0.1 


0.1440 4069 17 


0.2 


0.2276 3468 26 


0.5 


0.3213 4930 72 


1.0 


0.4511 8881 61 


2.0 


0.6192 0329 79 


4.0 


0.7179 6575 40 


6.0 


0.7455 6999 26 


7.0 


0.7582 2397 98 


8.0 



Table 4. Values of E^/m (ground state) for [ijm — 0.15 in 3+1. 



a 


Niewenhuis 
and Tjon [8b] 


Equal-time 
rcf. [8b] 


Gross eq. 
(with retard.) 


Present 
results 


Gross eq. 
ref. [8b] 


0.3 








1.999 536 




0.4 


1.99 






1.995 34 




0.5 


1.98 






1.986 30 




0.6 


1.96 


1.966 


1.969 


1.97176 


1.974 


0.7 


1.91 


1.941 


1.948 


1.950 81 


1.959 


0.8 


1.85 


1.907 


1.919 


1.92199 


1.938 


0.9 


1.77 


1.861 


1.880 


1.882 82 


1.910 


1.0 








1.828 47 




1.1 








1.746 64 




1.2 








1.562 48 




1.205 








1.534 20 




1.208 








1.503 66 




1.2085 








1.492 83 




1.2087* 








1.483 86 
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Table 5. Values of ^/m in 3+1 for fi/m = 0.15 for the first two excited states. 



a 


E 2 /m (n = 2,£ = 0) 


a 


E 2 /m (n = 2,£=l) 




(one node) 




(no nodes) 


1.0 


1.999 954 


1.4 


1.999 031 


1.25 


1.996 28 


1.5 


1.995 22 


1.5 


1.985 75 


2.0 


1.9518 


2.0 


1.9371 


2.5 


1.8461 


2.5 


1.822 7 


2.8 


1.686 2 


2.7 


1.725 7 


2.82 


1.6639 


2.8 


1.625 3 


2.84 


1.635 2 


2.82 


1.579 5 


2.85 


1.6161 


2.825 


1.556 8 


2.86 


1.589 2 


2.827* 


1.534 8 


2.866* 


1.5581 
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